Simple model for scanning tunneling spectroscopy of noble metal surfaces with 

adsorbed Kondo impurities 
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A simple model is introduced to describe conductance measurements between a scanning tunneling 
microscope (STM) tip and a noble metal surface with adsorbed transition metal atoms which display 
the Kondo efltect. The model assumes a realistic parameterization of the potential created by the 
surface and a d^,,2_^2 orbital for the description of the adsorbate. Fano lineshapes associated with 
the Kondo resonance are found to be sensitive to details of the adsorbate-substrate interaction. For 
instance, bringing the adsorbate closer to the surface leads to more asymmetric lineshapes while 
their dependence on the tip distance is weak. We find that it is important to use a realistic surface 
potential, to properly include the tunnelling matrix elements to the tip and to use substrate states 
which are orthogonal to the adsorbate and tip states. An application of our model to Co adsorbed 
on Cu explains the difi'erence in the lineshapes observed between Cu(lOO) and Cu(lll) surfaces. 

PACS numbers: 



I. INTRODUCTION 

The study of many-body phenomena in low dimen- 
sional systems is attracting a lot of attention. This has 
been motivated by recent advances in the construction of 
nanostructures and quantum dot devices. Scanning tun- 
neling microscopy has also opened the possibility of ana- 
lyzing many-body phenomena at surfaces. For instance, 
the Kondo effect has been detected in STM conductance 
measurements of noble metal surfaces with adsorbed 3d 
transition metal atomsiii2i2i^ Characteristic zero bias 
lineshapes are observed which are reminiscent of Fanc^ 
phenomena. 

Fano lineshapes have been observed in different situa- 
tions. They were first explored by FanoS. in his studies 
of autoionization of doubly excited He(2s2p) resonances 
lying in the continuum. In the case of adsorbed atoms on 
a metal surface an analogous situation is found as a lo- 
calized orbital is also coupled to a continuum of metallic 
states. As the adsorbate is magnetic, the Kondo effect 
can occur and, therefore, Fano lineshapes can be thought 
of as arising from the interference of the Kondo resonance 
with the continuum of metal states.^ It is worth noting 
that Fano phenomena appears in STM measurements al- 
though the tip overlaps much more strongly with the sur- 
face than with the adsorbate wavefunctions (as the 3d 
orbitals are very localized) so that conductance measure- 
ments reflect electronic properties of the metal surface 
modified by the presence of the magnetic atom. 

These Fano-type lineshapes differ from one adsor- 
bate/substrate system to another, as summarized in Ta- 
ble n although some experimental trends can be ex- 
tracted. For instance, the Fano parameter, q, is typically 
either zero or positive. As can be observed from the table 
the lineshape associated with g = is a symmetric dip 
close to zero bias. The dependence of lineshapes may be 
illustrated by comparing the asymmetric dip observed for 
Co on Cu(lOO) with the symmetric dip-like shape found 
for Co on Cu(lll). Similarly adsorbing Ti instead of Co 



TABLE I: Non-universality of lineshapes observed in conduc- 
tance measurements between a STM tip and a noble metal 
surface with adsorbed 3d transition metal atoms and Ce. De- 
pending on the adsorbate, the type of noble metal and/or 
the surface face, observed conductance lineshapes are difi^er- 
ent. The Fano parameter, q, which measures the degree of 
asymmetry of the lineshape associated with the Kondo reso- 
nance is typically g > 0. Ta' denotes the associated Kondo 
temperatures given in degrees Kelvin. 
Adsorbate/Surface Tk(K) Type of lineshape q Ref. 



Co/Cu(lll) 


54 


nearly symmetric dip 0.2 


[2,3] 


Co/Cu(100) 


88 


asymmetric dip 1.1 


f3] 


Co/Au(lll) 


75 


asymmetric dip 0.6 


fl] 


Ti/Au(lll) 


70 


asymmetric dip-f peak 1 


\8] 


Ti/Ag(100) 


40 


asymmetric dip-|-peak 1 


f4] 


Ce/Ag(lll) 


500 


symmetric dip 


f5] 


Co/Ag(lll) 


92 


symmetric dip 


m 



on Au(lll) leads to a strong variation of the lineshape. 
The situation becomes more complicated if we consider 
the " middle" elements of the 3d row which do not even 
show appreciable features in the conductancoSi down to 
T — 6 K. This raises questions about the occurrence 
of the Kondo effect at all for these specific adsorbates. 
Hence, experimental observations suggest that details as- 
sociated with the adsorbate-substrate interaction may be 
relevant. 

The Fano parameter, q, governing the shape of the 
conductance, G{uj), close to zero bias, is given by: 



AM 

where B{lu) reads: 

B(w) = 7r^Mkyk<5(w-ek), 



(1) 



(2) 



and, A{ll!), is the Kramers-Kronig transformation of 
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First attempts to model the substrate electronic struc- 
ture have been carried out by Plihal and GadzulsiS who 
have used a Jellium metal with a sharp step potential 
barrier at the surface. However, a half-filled symmet- 
ric density of states to describe the substrate (instead 
of parabolic) and momentum independent adsorbate- 
substrate hybridization matrix elements, Vk, were used. 
Under these strong assumptions dip-like conductance 
lineshapes were obtained. Tight-binding descriptions of 
the substrate have also been usediii They provide a qual- 
itative understanding of observations although their use 
is difficult to justify considering the sp free-like bands of 
noble metal surfaces. 



FIG. 1: Correspondence of shape of -B(aj) with the Fano pa- 
rameter. This schematic plot gives a qualitative understand- 
ing of how the shape of the conductance, G(tj), around the 
Fermi energy is related to B{ui). The Fermi energy refered 
to the vacuum level is denoted by the dashed vertical line 
e_F = —W, where W is the metal workfunction. — Vb denotes 
the bottom of the band. 

Bico): 

A{cu) = ^ r du^'^^. (3) 

In the above equations, Mk and Vk are the tip-substrate 
and adsorbate-substrate matrix elements, respectively, 
eic is the substrate band dispersion and ep the position 
of the Fermi level. 

In order to understand the relation between the en- 
ergy dependence of and the conductance lineshape 
we consider first in Fig. ^ two different situations at- 
tending to the degree of asymmetry of B{uo): (i) B(u)) is 
symmetric with respect to the Fermi energy. In this case, 
A{tF) — and the Fano parameter, g = 0. Hence, the 
resulting conductance lineshape is a symmetric dip. (ii) 
B{lo) is asymmetric with most of its weight situated at 
frequencies below the Fermi energy oj < ep- For this case, 
q > and conductance lineshapes are more asymmetric 
(q = f in the most asymmetric lineshape). Hence, ex- 
perimental observations requires that B{uj) is in between 
cases (i) and (ii). 

However, if we consider an sp-like band for the sub- 
strate and assume that hybridization matrix elements 
are independent of momentum, B{lj) would be parabolic 
(from Eq. lj2Jl), growing up to infinite energies. Due to 
the large asymmetry in B{uj), A{ep) would be very large 
and negative so that q ^ —oo leading to a peak (see 
Fig. ^ rather than the dip typically observed in experi- 
ments. Hence, it seems difficult to reconcile the observed 
conductance lineshapes with the shape of B(lo) from this 
simplified picture. The question that we address in the 
present work is whether the momentum dependence of 
hybridization matrix elements can change B{uj) leading 
to more symmetric conductance lineshapes. 



The above discussion shows the difficulty of describing 
experimental observations in a consistent way and points 
out the need of a more realistic model for describing the 
adsorbate-substrate-tip system. 

An issue which needs to be carefully addressed in defin- 
ing the relevant model is to know which of the states asso- 
ciated with the surface (either surface and/or bulk states) 
are more strongly coupled to the adsorbate. There is ex- 
perimental evidence suggesting that surface states play 
only a minor role in the metal-adsorbate interaction^. 
For instance, the amplitude of the conductance decays 
rapidly as, G cx with lateral displacement of the 

tip, instead of decaying as l/R\\ expected for surface 
states coupled to the impurity. Furthermore, this bulk- 
like behavior persists in conductance measurements of Co 
on Cu(lll), which is known to have a surface state at the 
Fermi energy^. Based on the above experimental obser- 
vations we will only consider bulk states in our model. 
Another experimental observation is that lineshapes de- 
pend only weakly with perpendicular tip-substrate dis- 
tance, Zt, suggesting that the direct interaction between 
the tip and the adsorbate is negligible 

In the present work we show how using an Anderson 
model in an appropriate orthogonalized basis, reasonably 
symmetric shapes of B{uj) are found which lead to con- 
ductance lineshapes in agreement with observations. We 
find that the momentum dependence of the hybridiza- 
tion matrix elements between the adsorbate Sd^^^-r^ 
orbital and the orthogonalized metal wavefunctions to- 
gether with the finite size of the tip wavefunction are 
responsible for this behavior. 

The present paper is organized as follows. In Section 
^we introduce an Anderson model in an orthogonalized 
basis to describe the adsorbate-substrate interaction to- 
gether with the relevant formulas needed for discussing 
STM conductance measurements. In Section fllll wc com- 
pute the parameters involved in the adsorbate- substrate 
interaction needed in the Anderson model. Section IIVI 
is devoted to describing the main results of our model. 
Finally, in Section we apply the model proposed to 
conductance measurements of Co atoms on noble metal 
surfaces. 
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II. THEORETICAL APPROACH 

In this section, details about the model used are given. 
An Anderson model on an orthogonal basis is consid- 
ered. The parameters of the model are computed taking 
into account the following considerations: (i) Instead of 
a step barrier we consider a Jones- Jennings-Jepsen (JJJ) 
potentialii to describe the metal surface wavefunctions, 
(ii) we neglect the direct coupling of the tip with the 
substrate d bands and with the 3d orbital of the adsor- 
bate due to the locaHzed nature of the d orbitals, (iii) 
the adsorbate is modelled by a single d-orbital, and (iv) 
the momentum dependence of the hybridization matrix 
elements is explicitly taken into account. 

Metallic states, |k >, coupled to a single d3z^-2 or- 
bital denoted by |c? > are considered. The Anderson 
model usually assumes that the cotinuum of metal states 
are orthogonal to the localized orbitals of the adsorbate. 
However, the basis set formed by the unperturbed metal, 
adsorbate and tip wavefunctions, {|k >,\d >,\t >}, is 
non-orthogonal (in general) and overcomplete. 

One way to take into account orthogonalization effects 
is to redefine the metallic states, k, as: 

|k >= |k > - < d|k > |d > - < t|k > |t > . (4) 

Considering that tip and adsorbate wavefunctions are or- 
thogonal: < t\d >= 0, as Id > is very localized, then the 
new metallic states satisfy 

<k|0>=O, (5) 

where |0 > can be either \t > or \d >. 

Our starting point is an Anderson model defined in 
this new orthogonal basis: {|k >, \d >, \t >}, from which 
associated one-electron parameters, eg., Vg. and Afj^ are 
obtained. An analogous procedure was previously used 
in the context of chemisorption of atoms and molecules 
on metal surfaces by Grimley^ 

It is worth mentioning that the new metallic states are 
non-orthogonal among each other: 

< k|k' >^ 0. (6) 

Hence, the orthogonalization condition between different 
|k > states is violated. However, this occurs at higher 
order in the overlap: 0{< kjd >^). 

A. Model 

From the above considerations, our model for the com- 
plete tip-substrate system (with the 3d transition metal 
atom) reads 

H = Hsubs + Htip-subs + Hfip, (7) 

where Hsubs describes the substrate with the adsorbed 
3d transition metal atom and Hup-subs describes the in- 
teraction of the tip with the substrate. Hup describes the 



tip which is assumed to have an unstructured density of 
states. 

The substrate with the adsorbed 3d atom may be mod- 
elled by a generalized Anderson model that explicitly in- 
cludes the orbital degeneracy of the 3d-orbital 

Haubs = ^ ^k'^^TOcr'^kmCT + ^ d]^^dm(y 
k,m,cr m,(y 

k,m,cj 

(8) 

Here is the energy level of an electron residing in the d 
orbital of the adsorbate, c- creates an electron with 

^ kmcr 

spin cr, momentum k and perpendicular projection of 
angular momentum, m, in the metal, d}^^ creates an 
electron in the state with perpendicular projection of an- 
gular momentum m in the adsorbate. ej. and are the 
metallic energies and the hybridization matrix elements 
between the substrate and the adsorbate, respectively. 
U is the Coulomb repulsion of two electrons in the 3d 
orbital of the transition metal atom. 

We will restrict the sum to m = as this is the or- 
bital which is more strongly hybridized to the metallic 
surfaccii^. In this case, Hamiltonian ((SJ reduces to the 
standard Anderson impurity model containing the spin 
degeneracy only. 

Finally, the tip-substrate interaction contribution to 
the Hamiltonian reads 

Hup-subs = E ^^k(4a*'^ + 
k,(T 

through the matrix elements, Afg.. Here, destroys an 
electron with spin cr in the tip. 



B. Hybridization matrix elements 

In the following we describe how hybridization matrix 
elements, Mg and Vg, are computed. For simplicity we 
will first focus on how Vg is computed, as Afg is computed 
in a similar way. 

Ignoring the electron-electron interaction we can re- 
express Hamiltonian ((SJ in first quantized form 

Hsubs ^T + Vd + VM, (10) 

where T is the kinetic energy of the system, Vd is the 
potential created by the adsorbate and Vm describes the 
surface potential. 

From Eqs. and IjlOl) . the matrix elements between 
the orthogonahzed metaUic states, k, and the adsorbate 
read: 

Vg-yk-^k<d|i4/M>, (11) 
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where, again, we have assumed: < t\d >= 0. 

The first term in Eq. is the hybridization matrix 
element with the unperturbed wavefunctions k; 

Vk=<k|^MM>, (12) 
and the second contains the overlap matrix element 



(13) 



The above orthogonalization procedure automatically 
selects the metal potential Vm in the hybridization ma- 
trix elements favouring the region close to or inside the 
metal in the integrations. This differs from hybridization 
matrix elements computed with the original wavefunc- 
tions as in that case integrations over the whole space 
are involved. 

Finally we note that orthogonalization effects enter the 
model through matrix elements only. Orthogonalization 
effects on the substrate band energies can be shown to 
be of higher order in the overlap. Hence, we will assume 
eg. = Ek in the rest of the paper. 

C. Computation of conductance 

Following Ref. ^3 and for the sake of clarity we derive 
the basic equations needed for the computation of the 
conductance through the STM. If we neglect any mod- 
ification of the substrate due to the presence of the tip 
(this is reasonable considering the fact that the tip is typ- 
ically at about 5 — 10 A above the metal surface), then 
the conductance measured by the STM reads ►ifi'ii'i^ 



Giu;) 



pup{T{uj) + ST{uj)), 



where 



(14) 



(15) 



is the conductance associated with the clean substrate 
(without the adsorbed 3d transition metal atom) and pup 
is the density of states of the tip. 

Modifications of the tip-surface coupling induced by 
the presence of the adsorbate are given by 



dr{u) = Im' 



k,k 



-Gdd[i^)- 



M*. Vt 

k' k' 



Ck' 



(16) 



In Eq. 1^3}, ptip is the density of states of the tip. 
Gdd{^) describes the electronic properties of the 3d ad- 
sorbate immersed in the metallic continuum including the 
many-body effects such as the Kondo effect induced by 
the on-site Coulomb interaction, [/, inside the localized 
id orbital, -q is an analytical continuation parameter. 

For convenience Eq. (|16|l is rewritten in the following 
way 

5T{lj) = lm{{A{uj) + iB{uj))Gdd{oj){A*{uj) + iB*{uj))} 

(17) 



where i?(w) is defined in Eq. (0) and A{u)) is the Kramers- 
Kronig transformation of B{ijj) given in Eq. Q with the 
matrix elements evaluated with the orthogonalized wave- 
functions, |k >. For the systems of interest here, Aiuj) 
and Biid) are real. Function B{ijj) embodies the informa- 
tion concerning the tip-substrate-adsorbate system as it 
depends on the tip-adsorbate separation, R, the position 
of the adsorbate with respect to the plane of ions, Zd, and 
the metal potential described by Vm through the matrix 
elements, Vg. and M^. In the present work we analyze 
how B[lij) and 5T[lij) depend on these parameters. 

From Eq. H17|l we notice that the conductance can, in 
principle, have any kind of shape as a result of the inter- 
ference of the adsorbate with the substrate continuum of 
states. In the following we will see how, in fact, Eq. H14|) 
reduces to the well known Fano expression. 



D. Fano formula for conductance 

For the sake of clarity we provide the Fano expression 
for the conductance which can be derived from Eqs. H14() - 
(|T7|l (details can be found in Ref. 0). First of all, the 
Green's function of the 3d orbital is assumed to be known 
and to have the simple form 



Gddiuj) 



Ad 



At 



Ak 



i2A 



uj — ep — £d — U ~ i2A 
(18) 



The three terms appearing in Eq. (|18|l correspond to the 
singly occupied atomic d level situated at e^, the Kondo 
peak situated at ck and the doubly occupied level, £d + U 
refered to the Fermi level, ep. A is defined as the half- 
width of the 3d impurity due to its hybridizaton with the 
metal surface: A{uj) = T^^]^\Vy.\'^S{uj — Ck) ■ Ak, Ad, and 
Ajj are spectral weights associated with the Kondo, the 
singly occupied and doubly occupied adsorbate, respec- 
tively. 

Greens function Gdd{oj) given by Eq. lfTH|l describes 
a magnetic impurity coupled to a metallic host in the 
Kondo regime. In the next section, we will give the pa- 
rameters relevant to the 3d impurities in noble metals 
used to model Gddi^^)- 

We introduce Eq. ^ in Eq. ^ and define e = 
For energies, uj < Tx and in the Kondo regime, T << 
Tk, the following Fano expression for the conductance is 
obtained: 



G(w) = C(w) 



(9 + ?) 



D{u), 



(19) 



where C{ui) and -D(w) have a weak dependence on uj and 
q denotes the Fano parameter given by Eq. Hence, 
the shape of the conductance observed, G{ijj), in the STM 
is governed by the value of the Fano parameter, q. Fig. 
n explains qualitatively the relation between B[u}) and 
conductance lineshapes. 



To conclude this section, we emphasize that Fano hne- 
shapes arising in the present model are a consequence 
solely of the interference between metal and adsorbate 
waves. No direct coupling between the tip and the 3d ad- 
sorbate is taken into account (which in any case should be 
very small) nor needed to explain Fano phenomena. How- 
ever, we will see below that the tip has, in fact, an im- 
portant influence on the final shape of B{lu) (and ST{uj)) 
through the matrix elements, Mj.. 



III. DISCUSSION OF PARAMETERS AND 
MODELLING OF THE SYSTEM 

From previous section it becomes clear that the shape 
of the conductance as a function of the bias measured 
by the STM depends on the parameters modelling the 
adsorbate-substrate-tip interaction entering B{lj). Here 
we give details of how these parameters are obtained. 

A. Surface potential and metal wavefunctions 

A realistic and at the same time simple description 
of Vm is attained by using a JJJ potential>i^ This po- 
tential interpolates between the image potential at long 
distances and the potential inside the bulk. It provides a 
realistic description of the surface barrier potential being 
particularly useful in interpreting LEED and photoemis- 
sion data of noble metal surfaces pi^^iSiii 

The JJJ potential contains three parameters that can 
be obtained from fitting the potential to density func- 
tional calculations or experimental data. The surface 
potential is given (in Rydberg energy units) by: 



Vm{Z) 



^1 (I - 

-Vo 



-A(Z-Zi„) 



) , z>z„ 



, Z < Zi, 



(20) 

Here Vq is the depth of the bulk potential and A controls 
the sharpness of the surface barrier potential. The pa- 
rameters A = 2Vb/A — 1 and /3 — Vo/A are obtained by 
imposing the condition of continuity of the potential at 
Z = Zim (see Ref. for more details). Throughout 
the paper we refer the tip, Zt, the adsorbate position, 
Zd, and the image plane position, Zim, to the last plane 
of ions situated at a (see Fig. ^ and all energies are re- 
ferred to the vacuum level. The sharpness of the potential 
barrier experienced by the electrons inside the metal is 
controlled by A. For increasing A, for instance, the po- 
tential becomes gradually sharper. As we will see below 
parameter A turns out to play an important role in our 
model. 

Metallic wavefunctions are obtained considering box 
quantization. We define a large box of side 2a, which 
describes the substrate enclosed by a larger box of side 
2L, where the relative sizes of the boxes satisfy: L >> a 
(see Fig. In practice it is sufficient to take a/L = 0.3. 
The perpendicular dependence of the wavefunction is ob- 
tained from the integration of Schrodingers equation in 
the presence of the JJJ potential starting at distances of 




FIG. 2: Schematic representation of the potential used to 
model the noble metal surface. The sizes of the large and 
small boxes used for quantization of the wavefunctions are 
given by 2L and 2a, respectively. The hatched vertical rect- 
angles denote walls of infinite potential. Zd and Zt denote 
the positions of the adsorbate and the tip with respect to the 
last plane of ions. Vo is the height of the surface barrier and 
W the metal workfunction. The Fermi energy is referred to 
the vacuum level and is denoted by the horizontal dotted line: 
ep = -W 



about 25 A outside the surface towards the metal. Impos- 
ing periodic boundary conditions in the parallel direction, 
the crystal wavefunction reads: 



*(r) 



„«k||r 



n 



(21) 



where, ^I'fe^ (z), is the perpendicular part of the full wave- 
function. Sates above the vacuum are normalized to the 
volume n = while states below it are normalized to 

B. Tip and adsorbate wavefunctions 
The adsorbate wavefunction is described using a 
Slater-type function: 



Mr) = Ny-'e 



n — l —rjir 



Yin 



(22) 



where n is the main quantum number, / the angular 
momentum and m the z-component of the angular mo- 
mentum. The normalization constant of orbital i is 
Ni = [(2n)!]-i/2(2,yi)("+i/2) and Yim{e,(j)) is a conven- 
tional spherical harmonic. For the 3d transition metal 
atom we use n = 3, ^ = 2, m = with an exponential 
decay rjd = 2.^0^^ , where oq is the Bohr radius. 

Following Tersoff and Hamann,^^ the tip wavefunction, 
\t >, is described by a simple s-type orbital [n = 1), with 
an exponential decay fixed by the workfunction of the 

metal surface: rjt = \J 2mf,W/h? , where vrif. is the elec- 
tron mass and W is the workfunction of the metal surface 
which we take to be W = 4.95 eV in our calculations. 
This choice for rjt relies on the fact that, typically, tung- 
sten tips are dipped in the metal surface before measuring 
the conductance between the STM tip and the surfacei^ 
Hence, it is reasonable to describe the wavefunction tails 
sticking out from the tip with the same exponential decay 
as the ones describing the metal surface. 
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C. Computation of matrix elements 

Once we have obtained the wavefunctions of the sub- 
strate we may proceed with the computation of the 
hybridization and overlap matrix elements defined in 
Eq. Hll() . For simplicity we discuss the adsorbate- 
substrate hybridization matrix elements. From Eq. 1)2 l|l 
hybridization and overlap matrix elements expressed in 
cylindrical coordinates, read: 

/•OO POO 

= 2tt dr\\rn / dz^-fc, (z) Jo(fc||r||)VA,/(z)0d(|r - 

JQ Jo 

/>00 /"OO 

S'k = 27r / dr||r|| / dz*fe^(z)Jo(fc||r||)(/)<j(|r - rd|) 
Jq Jq 

where Jo(fc||r||) is the cylindrical Bessel function of zeroth 
order, Vm, is the surface potential and (^d(|r — i"d|) is the 
adsorbate wavefunction. The computation of the tip- 
substrate hybridization matrix elements, Mk, proceeds 
along similar lines replacing (j)d{\r — r^l) by 0t(|r — R|). 

The behavior of tip matrix elements deserves careful 
attention. Due to the slow decay of the tip wave function, 
most of the contribution to comes from the region 
close to the metal surface for not too large. Matrix 
elements Afj.(R) can then be approximately factorized 
asiS: 

M£(R)«/(Zt)e^''i|R|iAfi^. (24) 

where Mj. is independent of the tip position. This re- 
mains true for kz ^ 1.6 — 1.7 A~^, which is close to 
ko = 1.85 A-i (for Vo = 13.05 eV). The function f{Zt) 
is related to the value of the tip wave function in the sur- 
face region and it decays rapidlyi^ with Zt . Mj^ is found 
to have a weak dependence on kz in contrast to matrix 
elements obtained from the original substrate wavefunc- 
tions as can be observed in Fig. 13 Indeed, using Va/ 
in the matrix elements as imposed by the orthogonaliza- 
tion cuts off part of the exponential dependence of metal 
wavefunctions far from the surface leading to a weaker 
dependence on kz. This is crucial to obtain reasonable 
shapes of B{uj). 

In order to have a qualitative understanding of the 
behavior of adsorbate-substrate hybridization matrix el- 
ements, Vk, we focus on electrons which have no parallel 
momentum, fc|| =0. Matrix elements gradually increase 
with increasing kz due to increasing overlap of the metal 
states with the adsorbate. The behavior of matrix ele- 
ments at large kz can be partially understood by analyz- 
ing Eq. (|23|l with ^'^^(z) nearly constant, (which is the 
case for the metal wavefunction outside the surface). For 
a d^z^-r^ adsorbate the associated matrix elements, 
would tend to zero due to cancellations associated with 
the lobes of the 3d-orbital. Indeed, we find that while 
Vj^ computed for a d-orbital presents a downturn at in- 
termediate kz ~ 1.4 A~^, computed with an s-type 
orbital increases rapidly up to kz « feg. Orthogonaliza- 
tion effects are also important as can be noticed from the 
fact that Vk ~< d.\VM\d > 5k as kz ko. Hence, we 



0.002 



0.000 




0.0 0.5 1.0 1.5 

K 



FIG. 3: Effect of orthogonalization on tip-substrate hybridiza- 
tion matrix elements. While hybridization matrix elements 
computed between the original substrate wavefunctions and 
the tip depend strongly on perpendicular momentum, fcz, this 
dependence is weaker for matrix elements containing orthogo- 
nalization corrections. The JJJ potential used in this calcula- 
tion is parameterized using A = 1.7 A~^, Vq — 13.05 eV, and 
Zim = 1.1 A. The tip is located at Zt — 5 K above the surface 
plane of ions on top of the adsorbate (i?| | = 0) . kz is given in 
with the vacuum level corresponding to fco — 1.85 A~^. 



find that matrix elements computed using orthogonal- 
ized wavefunctions, l/j., initially increase with kz, reach a 
maximum and then they are gradually suppressed regard- 
less of the angular dependence of the adsorbate orbital. 
This dependence of Vg. is reflected in the frequency de- 
pendence of B{uj). The above discussion illustrates the 
importance of appropriately including orthogonalization 
effects in the computation of hybridization matrix ele- 
ments. 

The behavior of matrix elements at very large per- 
pendicular momentum needs to be carefully examined. 
Metallic waves with large momentum (corresponding to 
energies far above the vacuum level), contain rapid os- 
cillations which lead to cancellations in the integrand of 
Eq. (|23|l . A high energy cut-off may be simply estimated 
from: kf^'^°^^ ~ rji. This gives lo'^^^°^^ « 70 eV for 

the adsorbate and uj^"*""^^ ~ 5 eV for the tip. Hence, 
hybridization matrix elements and in turn B{uj) will be 
non-zero in the energy range — Vq < w < w™*"^-^. It is 
remarkable that the tip cuts off all high energy contribu- 
tions having an important influence on the final shape of 
B{ijj) and in turn on conductance lineshapes. This point 
seems to have been overlooked before. 
D. Computation of B{uj) 

Once we have the hybridization matrix elements, it is 
straightforward, following Ref. to obtain the energy 
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dependence of B{uj) given in Eq. After performing 
integration in parallel momentum, B[u) can be written 
as 

B{u;)^-^-^{u^ + Vo) dxVj^^^;,M~,^^~,^, (25) 



where k\\ = y 2mf./h?{uo + Vq)\/1 — and = 

\J 2mi,/'li' [lo + Vb)a; and is the volume. Note the factor 
il appearing in front of the integral cancels out the vol- 
ume coming from the normalization of the surface wave- 
functions coming in the matrix elements. 

E. Parametrization of adsorbate Greens function 

In order to compute the modification induced by the 
adsorbate on conductance lineshapes, 5T{uj)^ knowledge 
of Gdd{'-^) is needed. We take parameters used by 
Ujsaghy et. ali to model Co on Au: = 3 meV, 
Ed = -0.84 eV, U = 2.8 eV, Ak = ttTk/A, Ad = 0.45 
and Au = 1 — Ad~ Ak with the Kondo temperature fixed 
at Tk — 5 meV. As we are particularly interested in what 
is happening close to the Fermi energy (at energy scales 
of the order of Tk) the relevant parameters to our calcu- 
lations are Tk and ck- Hence, these are input parame- 
ters which could be obtained from other sources such as 
experimental data or ab initio density functional calcula- 
tions. It is not the purpose of the present work to have a 
first-principles determination of Kondo temperatures for 
different adsorbate/substrate systems but rather to have 
a qualitative understanding of conductance lineshapes. 



IV. DISCUSSION OF RESULTS 

In this section we show our main results and discuss 
their relevance to experiments. A positive Fano param- 
eter g > is typically found in our calculations for pa- 
rameters charaterizing the noble metal surface, the ad- 
sorbate and their mutual interaction. We explore how 
lineshapes depend on the parameters characterizing the 
atom-surface interaction. These parameters are: the po- 
sition of the tip R = (R||,Zf), the adsorbate-substrate 
distance, Zd, and the shape of the surface potential con- 
trolled by A. As there is cylindrical symmetry in our 
model, our results only depend on the absolute magni- 
tude of the tip-adsorbate lateral displacement — |R|||. 
We do not find much dependence of our results with the 
workfunction or the height of the surface barrier so we fix 
Vq = 13.05 eV and the workfunction W = 4.95 eV, which 
are typical values for noble metal surfaces. We will first 
discuss the dependence of lineshapes with the tip fixed 
on top of the adsorbate, = at a distance of about 
Zt = 5 A. At the end of the section we will study how 
lineshapes depend on and Zt- 



A. Dependence of Fano lineshapes with 
adsorbate-surface distance 

Fixing the tip at =0 (on top of the adsorbate) 
and Zt = 5 A we analyze how lineshapes depend on the 
adsorbtion distance, Zd- In Fig. 21 the dependence of 
B{uj) and ST{uj) is shown for typical parameters describ- 
ing noble metal surfaces: A ~ 2.08 A^^, Vq ~ 13.05 eV 
at typical adsorbate-surface distances: Zd ~ 1.5 — 2 A. 
Adsorbtion distances for different surfaces are obtained 
assuming hard spheres for the atoms. At large adsorb- 
tion distances, Z^ « 2 A, B{uj) is found to be rather 
symmetric with respect to the Fermi energy and conse- 
quently (see Fig. 2) the lineshape is a nearly symmetric 
dip. As the adsorbate is moved closer to the surface 
{Zd decreasing) the frequency dependence of B{uj) varies 
and consequently the lineshapes too. An enhancement of 
weight below the Fermi energy as well as a sharper drop 
above it occurs. This leads to a variation of lineshapes 
from g « at large distances towards g > at short 
adsorbate-surface distances. 

To understand these results, we notice that substrate 
states with low energy have a strong exponential decay 
outside the surface. As the adsorbate is brought closer 
to the surface, the overlap of the adsorbate orbital with 
these states increases. This leads to an increase of B{lu) 
at small tu. States with a larger energy extend farther 
outside the surface and have a substantial overlap with 
the different lobes of the adsorbate 3c?-orbital. This tends 
to lead to a cancellation of the different contributions to 
Vk. This tendency increases as the adsorbate is brought 
closer to the surface, explaining the reduction of B{uj) at 
larger uj. 

For comparison we also show in Fig. [Slthe dependence 
of lineshapes with the adsorbate-surface distance for a 
Jellium with a step potential. Lineshapes are found to be 
somewhat asymmetric as a large part of the weight in sit- 
uated at low energies in B{lu) (see also Fig. A weaker 
dependence of B{uj) with the adsorbate-surface distance 
is obtained as compared to the JJJ potential. The fact 
that we obtain reasonable lineshapes with a simple model 
such as a Jellium with a step potential can be attributed 
to orthogonalization effects. Tip-substrate matrix ele- 
ments computed without orthogonalization corrections 
would increase rapidly with kz as they involve integra- 
tions over the whole space (instead of the region close 
to the surface) leading to rather asymmetric shapes of 
B{u;). 

B. Dependence of lineshapes with the shape of the 
potential barrier 

There are two parameters that can be tuned to change 
the shape of the surface potential. One is the image 
plane position, Zim, and the other is the sharpness of 
the surface potential given by A. We fix = 1.1 A and 
study the dependence of -B(a;) and lineshapes Sr{uj), with 
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FIG. 4: Dependence of lineshapes on the adsorbate-metal ad- 
sorption distance. As the adsorbate is brought closer to the 
surface the hneshapes become gradually more asymmetric. 
The lower pannel displays the function -B(a;) and in the up- 
per one is the corresponding modification in the tip-surface 
interaction induced by the adsorbate: Sr{u)). The JJJ po- 
tential is parametrized using A = 2.08 and Zim = 1.1 
A. The tip is located a,t Zt — 5 K above the surface plane of 
ions. All distances in the figure are given in A and refered 
to the last plane of ions. The vertical dotted line denotes the 
Fermi level position: ep ~ —W — —4.95 eV. The height of 
the surface barrier is Vb = 13.05 eV. 



the sharpness of the potential barrier. This dependence 
is shown in Fig. |H| We find that an increase in A leads to 
a more asymmetric B{lu) and the associated conductance 
lineshape, 5T{u}) with more positive values of the Fano 
parameter, q. 

Increasing A makes the metal potential, Va/, sharper so 
that the exponential tails of the metal wavefunctions out- 
side the surface are shifted closer to the adsorbate leading 
to a stronger metal-adsorbate overlap. This is analogous 
to bringing the adsorbate closer to the surface for a fixed 
A as discussed previously. Hence, an increase in A leads 
to an increase in the weight of i?(ti;) at low energies and 
a sharper drop at energies above the Fermi level. The 
stronger dependence found for B{uj) with A than with 
the adsorbate-substrate distance can be ascribed to the 
exponential dependence of Vm with A. 

Summarizing, we find that lineshapes are strongly de- 
pendent on the degree of sharpness of the metal poten- 
tial. Making the surface potential sharper leads to more 
asymmetric lineshapes at typical atom-substrate adsorp- 
tion distances. 



FIG. 5: Dependence of lineshapes on the adsorbate-metal ad- 
sorption distance for a Jellium model with a step potential. 
In contrast to the JJJ model potential (see Fig. |^ there is a 
weak dependence of lineshapes with the adsorbate-substrate 
distance. The tip is located at Zt — 5 K above the sur- 
face plane of ions. All distances in the figure are given in 
A. The vertical dotted line denotes the Fermi level position: 
ep = —W — —4.95 eV. The height of the surface barrier is 
Vo = 13.05 eV. 
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FIG. 6: Dependence of lineshapes with the sharpness of the 
surface potential. As A increases B{uj) and lineshapes become 
more asymmetric. The lower pannel shows B(lu) for different 
A and the upper one the corresponding lineshapes. The sur- 
face potential is modelled for different A and Zim = 1.1 A. The 
position of the tip and adsorbate are Zt — 5 K and Za — 1.6 
A, respectively. All distances are given in units of A and A is 
in units of A~^. The vertical dotted line denotes the position 
of the Fermi level: ep = -W = —4.95 eV. 
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C. Dependence of conductance lineshapes with tip 
position 

We finally discuss the dependence of conductance line- 
shapes when the tip is moved away from the adsorbate 
either laterally 7^ 0) or perpendicularly. The tip 
position comes into our calculations through matrix ele- 
ments, Mj^. A weak dependence of lineshapes with per- 
pendicular tip-substrate distance, Zt is found. Varying 
the distance from 5 to 7 A leads to a small change in the 
lincshape becoming slightly more asymmetric (see Fig. 
[Tj). The weak dependence of the line shape on the tip- 
surface distance can be understood from Eq. (|24|l . This 
equation predicts a strong reduction of the amplitude 
of the conductance with vertical displacement of the tip 
^ f{Zt)'^, but no change in its shape. The weak depen- 
dence found is consistent with experimental observations 
which also indicate that the direct tip-adsorbate coupling 
is very small (as expected from the localized nature of the 
d-orbital of the adsorbate). 

We can also analyze how lineshapes depend on the lat- 
eral position of the tip, An example is shown in Fig. 
[3 which a rapid decay of the amplitude of the lineshapes 
is found. Lineshapes are found to depend slightly on 
at tip-surface distances of about Zt = 7 A which is the 
typical experimental tip position. Indeed, crystal wave- 
functions with A:|| = have the largest amplitude outside 
the surface representing electrons which have the largest 
probability of being detected by the tip (see Ref. for 
discussion). Hence, the oscillatory behavior appearing in 
Eq. H24|) is suppressed as only the states with small 
contribute to the momentum sums. 

The rapid decay of the amplitude with obtained 
from our calculations is consistent with the disappear- 
ance of the features in the conductance at distances of 
about ~ 5 — 10 A. However, the small variation of 
the shape of the conductance with found is only par- 
tially consistent with experiments as for some systems 
such as Co on Cu(lOO) or Co on Au(lll) lineshapes be- 
come gradually more symmetric {q — > 0) at larger R^. It 
is worth mentioning that in the case of Co on Au(lll), 
the lateral dependence of lineshapes depends not only on 
the modulus of R|| but also on the direction in which 
the tip is moved^ from the adsorbate. This cannot be 
described in the present framework as the surface corru- 
gation is not included in our model. 

The finite size of the tip leads to a suppression of the 
tip-substrate matrix elements at large energies due to 
rapid oscillations of the metal wavefunctions providing a 
high energy cut-off. This cut-off turns out to be crucial 
to obtain reasonable lineshapes. Had we used a point-like 
tip the amplitude of B{uj) at energies below the vacuum 
level would be exponentially suppressed while above it 
large oscillations would occur. This is because wavefunc- 
tions are evaluated far from the surface (at the tip po- 
sition). Hence, for a point-like tip the Fano parameter 
would be very large {q >> 1) and a peak rather than a 
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FIG. 7: Dependence of lineshapes on the position of the tip. 
The upper panel shows the small perpendicular dependence 
of lineshapes if the tip is moved outwards from the surface 
by 2 A. The surface potential is parametrized by A — 2.1 
A~^ and the position of the image plane Zim — 1.1 A. The 
adsorbate is at Zd = 1.6 A. All distances in the plot are given 
in A and the Fermi level is situated at: ep = —4.95 eV. The 
conductance lineshape at = 7 A is multiplied by 40 for 
comparison with Zt = 5 K. 

dip would appear in conductance measurements. 

V. AN APPLICATION TO CONDUCTANCE 
MEASUREMENTS OF COBALT ON COPPER 
SURFACES 

In the previous sections a detailed description of the 
method used has been provided. In the following we pro- 
vide the reader with an application to conductance mea- 
surements of Cu surfaces with Co atoms deposited on 
them. 

Copper surfaces are described using the parameters 
quoted in Ref. jl3| for the analysis of surface states in 
noble metal surfaces. For Cu surfaces we take A = 2.2 
A~^ and the adsorbate -substrate distance is estimated 
assuming hard spheres for the atomsiSS As the (111) sur- 
face is more closely packed than the (100) surface the 
Co atom adsorbs at larger distances in the former than 
in the latter surface. This leads to Zd — 1.5 A and 
Zd = 1.8 A for the (100) and (111) surfaces, respec- 
tively. At the temperatures used in the experimental 
studies of the conductance?^ Co sits outside the Cu(lOO) 
surface and it is not incorporated in the surface, as hap- 
pens at higher temperaturesSl. Ab initio calculationsSS 
then give Zd = 1.5 A, in agreement with our simple esti- 
mate. We find that difference in the adsorbtion distance 
is enough to produce changes in the lineshapes as can 
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FIG. 8: Lineshapes associated with conductance measure- 
ments of Co on Cu(lOO) and Cu(lll). Parameters taken from 
Ref. [l3 are used to describe the surface potential: A = 2.2 
and Zim = 1.15 A. The adsorbtion distance is larger for 
Cu(lll) than Cu(lOO) by about 0.3 A. This leads to a more 
asymmetric lineshape in the latter than in the former case in 
agreement with experimental findings.'^ All distances in the 
figure are given in A and A is given in A~^. 



be observed in Fig. |S1 From the dependence of line- 
shapes with Zd discussed above we find that for Co on 
Cu(lOO) surfaces lineshapes would be more asymmetric 
than for Co on Cu(lll). This is consistent with experi- 
ments performed for these two different surfaces^*^ Our 
result may be of more general validity and may apply to 
conductance measurements on other metal surfaces such 
as Ag or Au. However, care should be taken as our model 
contains the adsorbate-substrate distance only: corruga- 
tion effects may also be important. More experiments 
that systematically analyze the dependence of lineshapes 
for Co on different crystal faces are needed to corroborate 
our findings. 



For typical values of the parameters characterizing the 
adsorbate-substrate interaction, we find lineshapes with 
Fano parameters, g > 0, in agreement with experimental 
trends. 

We have introduced a function B{uj) which describes 
the adsorbate-substrate interaction. The description of 
experimental observations requires a B{uj) which is fairly 
symmetric with respect to the Fermi energy. Since the 
underlying density of states is approximately parabolic 
and the matrix elements tend to further reduce B{ll}) at 
small uj, it is interesting that the model, nevertheless, 
gives reasonable shapes for B{uj). We find that the tip 
matrix elements play an important role. The tip wave 
function is rather extended, and matrix elements to con- 
tinuum states at high energies are very small, due to the 
rapid oscillations of the continuum states over the spa- 
tial range of the tip wave function. This helps to make 
B{ll!) more symmetric around uj = ep. We also find that 
it is essential to use substrate states which have been or- 
thogonalized to the adsorbate and the tip. Even when 
the realistic JJJ potential is used, lineshapes are found 
to strongly depend on the sharpness of the surface po- 
tential controlled by A becoming more asymmetric with 
increasing A (i. e. for sharper potentials). To a lower 
degree conductance lineshapes are also found to depend 
on the adsorbtion distance: the closer to the surface the 
adsorbate is the more asymmetric the lineshapes become. 
These tendencies can be explained from the enhancement 
in the coupling of the metal wavefunctions and the ad- 
sorbate at low energies together with the suppression at 
energies above the Fermi level induced by the angular 
dependence of the c?3j,2_^2 orbital and the orthogonaliza- 
tion. 

More closely packed surfaces such as Cu(lll) are ex- 
pected to lead to more symmetric lineshapes than the 
more open Cu(lOO) faces as Co adsorbs further out in 
the former than in the latter. We find good agreement 
with experiments comparing lineshapes of Co on Cu(lOO) 
with Co on Cu(lll) surfaces. In agreement with exper- 
iments we also find that the tip-substrate distance does 
not strongly influence the shape of the conductance. 



VI. CONCLUSIONS 
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